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The low energy effective field model for the multilayer graphene (at ABC stacking) in external Electric 
field is considered. The Schwinger pair creation rate and the vacuum persistence probability are calculated 
using the semi - classical approach. 
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1. INTRODUCTION 



The Schwinger mechanism of electron - positron pair 
creation in electric field was first investigated in [T] . The 
rate of pair production depends on electric field and is 
so small that is not observable for the experimentally al- 
lowed values of the electric field. In the condensed mat- 
ter systems the situation may be different. For example, 
in graphene [31 [5] the large value of the effective coupling 
constant opens the possibility for the pair creation pro- 
cess to be observed [H [S] . In graphene monolayer the 
pair creation rate may be calculated using the approach 
■of [I]. However, a different approaches were also used 
(see, for example, [U E] and references therein). The 
approach described in [B] was also applied to the bilayer 
graphene [HIS]- 

In the present paper we calculate the pair produc- 
tion rate in multilayer graphene. We consider the sim- 
plest case of ABC - stacking described by the two - 
band pseudospin Hamiltonian with the chirality index 
J equal to the number of layers [TU]. In our calcula- 
tions we rely on the method developed in [TTl |3] and 
in [12]. A similar approach was also applied to IIe-3 
[7]. Within this approach, which was originally applied 
to monolayer graphene, we develop the semi - classical 
approximation (for the alternative ways to apply semi - 
classical technique to the fermionic models see [T31ll4j). 
This approximation gives results identical to the results 
obtained via the exact solution of the Schrodinger equa- 
tion for the case of monolayer. Our results obtained in 
the multilayer graphene are checked with the more tra- 
ditional semi - classical approach described in [BJ [HI [H] ■ 
The paper is organized as follows. In Section 2 we 
describe the one - particle Schrodinger equation that 
appears in the given problem. In Section 3 we intro- 
duce appropriate boundary conditions. In Section 4 
the semiclassical approximation for the one - particle 
Schrodinger equation is introduced. In Section 5 we 



check the results obtained in Section 4 via the semiclas- 
sical approach described in [6j HI [9]. In Section 6 we 
compare our results with the exact ones for the case 
of graphene monolayer. In Section 7 we calculate the 
vacuum persistence probability and the pair produc- 
tion rate for the field - theoretic model of multilayer 
graphene. In Section 8 we end with the conclusions. 



2. ONE - PARTICLE SCHRODINGER 
EQUATION 

First, let us consider the one - particle problem. We 
deal with the two - component spinors placed in the 
external Electric field directed along the x - axis. We 
consider the external Electro-magnetic potential in the 
form: = Et. The one - particle Hamiltonian in a 
subsequent parametrization has the form |101 [5] 
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(1) 

Here w is a constant that is equal to Fermi velocity 
for the case of monolayer. J is the number of layers. 
Schrodinger equation has the usual form 



(2) 



Its solution is 
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^it) = Pexp (-1 j H{t)dt\^{to) = Umito) (3) 

J to 

Here the path - ordered exponent is used. Operator U{t) 
is unitary by construction. Later on we imply periodic 
boundary conditions in space - coordinates. That's why 
^ can be decomposed into the sum over the quantized 2 
- momenta: ^{t,x)) = Ep,,p„ e'P-''+'Pyy^pp{t). For V'(t) 
we have: 
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<p{t) = Pexp (~*'^^ H\p,t\dtjtljp{to 



(4) 
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3. BOUNDARY CONDITIONS 

In our semiclassical consideration we imply that J 
is odd. However, analytical continuation will allow us 
to obtain final results for even values of J as well. It 
is implied that at t < to Electric field is absent and we 
have 



i;p{t) = R\p,t]+e 



-iv\P\-'(T3t 




(5) 



where P = Px + Eto + ipy Boundary conditions at to 
must correspond to the negative energy levels occupied: 



(6) 



It is supposed that at t > to + T Electric field is 
switched off again. Then at t = to + T we have 

'^-^^[V'^-QV] (7) 

the value |?7+p is the probability that the electron - hole 
pair has been created, while j??- p is the probability that 
the negative energy level remains occupied. 

For the semi - classical consideration it is useful to 
consider to = —T/2 — )• — oo. The essence of semi - 
classical methodology is the consideration of large fre- 
quencies that are in this case Px + Et. That's why we 
require at t = to = —T/2 



= V-'l - -02 



(8) 



Att = +T/2 we denote 
1 



V2 
1 
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(9) 



where again the value |?7+P is the probability that the 
electron - hole pair has been created. This consideration 
implies -ET/2+px < and ET/2 +px>0. 

For ET/2 < px (px < -ET/2) boundary conditions 
att = to = -T/2 are: 



= Vl ± ^-2 

1 = ;^[^2TV'i]- 



(10) 



Here the upper sign is for ET/2 < Px while the lower 
one is for -ET/2 > px- A.tt = -\-T/2 we expect 



??+ = [V'l ± V-'2] 

n- = -^[02 T V'l], 



(11) 



where again the value |?7+p is the probability that the 
electron - hole pair has been created. 



4. SEMICLASSICAL CONSIDERATION 

Let us now introduce the notations: 

r = (^)*fc + £.), n = (|)*^, 

Q = {T + iJiy (12) 



Then 



il)p{t) = Pexp 



e* 
e 



dT 



)Mto) (13) 



The corresponding system of equations at t > to is: 



itp2 = Qtpi 
For Vi,2 we have: 

il)2 = iip'i/w^ , u = T — iIl 



(14) 



l-jduj^iJi + (1 + — )^i = 



2m, 

u 



(15) 



Wc introduce new variable z = u'^~^^/{J + 1). The 
resulting equation is 

[^^]'^^ + (1 + (jTi)otWtT)^'^^ = ° 

We represent tpi = e** and obtain the equation for 
s that is considered iteratively. In the first approxima- 
tion we neglect the derivatives of s higher than the first 
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derivative. In order to calculate the second approxi- 
mation we substitute the first approximation to the ex- 
pression for s" etc. When only the first and the second 
terms are kept, the wave functions are given by 



+C2 



e '^0 
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The considered approximation is valid if the second 
approximation is smaller than the first one. This leads 
to the condition 
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Next, we use the fact that the given semi - classical 
approximation gives the solution of Eq. (llSp not only 
for large real values of r but for the complex values of 
T with large |r|. That's why analytical continuation of 
the solution at t — — oo gives the solution at r ^ +oo. 
(The continuation is performed along the line placed at 
|r| -5> c».) 

At ET/2 > px > —ET/2 boundary conditions give 
C2 = 0. The probability that the pair is created is |?7+P, 
where 

^^^\ni'+^BiJ/2+i.i/2)^ z^-(r/n)2 (19) 

Here we consider the contour placed at infinity with the 
orientation such that \ri^\ remains less than unity. 
We have 
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(20) 



Written in this form our result can be continued an- 
alytically to even values of J . At the same time known 
results for J = 1,2 are reproduced [H E]- 

At ETI2 < px {px < -ET/2) boundary conditions 
give ci = (c2 =0). In both cases semicassical approx- 
imation gives |?7_| = 1 that means that the electron - 
hole pair is not created. 



Below we check the obtained above value of the prob- 
ability that the electron - hole pair is created with the 
given values of momenta (px^Py). We do this in two 
ways: via the application of the semi - classical approx- 
imation in its more classical form and via the consider- 
ation of the exact solution of the Schrodinger equation 
(at J = 1). 

5. MORE CLASSICAL SEMI-CLASSICS 

The problem of pair creation can be considered 
in the gauge Aq = —Ex. Then we have stationary 
Schrodinger equation H'ii = e^P with 




' [Px + iPy j 

, X 1/(J+1) 

We proceed with the rescaling z — l-^j x, and 



1/(^+1) 



e. Then 



{z - uj)ipi + {-id, - ffl)'V2 = 
(z - a;)V'2 + [-id., + jn)^Vi = 



(22) 



The first order semi - classical approximation for "^x^i 
gives 



(23) 



Integration over the classically forbidden region —112 _|_ 
(z — uj)"^! < gives the pair production probability: 

a^2.JB{l,l) (24) 
This expression coincides with the one derived above. 

6. EXACT SOLUTION AT J = 1 

Let us introduce notations ip+ = ipi — 'ip2 and 
■0_ = ipi — ip2- Then 

1 T 

V-^ + (i + n2 +t2)V+ = (25) 
We change variables r = -^e^'^^^z. Then 

V'- = ^-4t [d,-z/2]^j+ 



n 



+ (1/2 + - 1] - zyA)ij+ - (26) 
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The solution is 



^+ = 6ii?_,n2/2(y2e-3-/V) +62AnV2-i(V2e3-^/V) 



{-&inV2i?-ffl2/2_i(\/2e"3-/V) 



n 

-62AnV2(V2e3-/V)} 



(27) 



The consideration of usual boundary conditions 
leads to rather complicated algebra. So we come to the 
semiclassical boundary conditions (tq — > +00): 

V2 = 6ii?_,nV2(^e"^^^o) + 52AnV2-i(V2e-^"/Vo) 
= -6inV2i?_,n2/2-i(\/2e^"/Vo) 

+62A:nV2(V2e-"/Vo) (28) 

Using asymptotic expansion for Weber function we 
come to 62 = and 



(29) 



is the probability that the electron - hole pair is created. 



7. FIELD - THEORETICAL CONSIDERATION 

The fact that the particles do not interact with each 
other allows to reduce the field - theoretical problem 
to the quantum mechanical one. Namely, we arrive at 
the following pattern. Modes for different values of mo- 
menta propagate independently. At i < states 
with negative values of energy are occupied while all 
states with positive values of energy are vacant. Their 
evolution in time is governed by the one - particle 
Schrodinger equation. At t = to -\- T the wave func- 
tion already has the nonzero component corresponding 
to positive energy. Its squared absolute value is the 
probability that the electron - hole pair is created. 

Let us calculate the probability that vacuum remains 
vacuum (vacuum persistence probability). Accord- 
ing to the above presented calculation this probability 
is 



n.,4i-e-("-"^"""')"")' 



,-2IinS 



(30) 

Here S is the effective action, the factors gs — 2 and 
Qy = 2 are spin and valley degeneracies. The product is 
over the momenta that satisfy 

ET/2 >p^> ~ET/2 (31) 

We have [L is the linear size of the graphene sheet): 
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~ ^'^^ 2( J + l)7r2(a)i/(./+i) ^^^^ 

Here a is given by Eq. ([20]). According to a 
different quantity is considered as the pair production 
rate: 
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(33) 



The form of the functional dependence of F on E 
coincides with that of mentioned in [9]. 



8. CONCLUSIONS 

In the present paper we calculate the pair produc- 
tion rate and the vacuum persistence probability for the 
multilayer graphene. We develop the semiclassical tech- 
nique within the approach used earlier in monolayer 
graphene. Our method reproduces known results for 
monolayer and bilayer graphene. Following jlll we con- 
sider the single pair creation rate F and uj 



(where P^ is the vacuum persistence probability) as dif- 
ferent quantities. The possibility to consider w as a 
production rate of multiple states remains open and re- 
quires an additional investigation. 
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